1 



ON THE NATURE OF THE UNIDENTIFIED MEV GAMMA-RAY SOURCE GRO J1411-64 



Gustavo E. Romero 1 , Mariana Orellana 1 , and Diego F. Torres 2 

1 IAR, C.C. 5, 1894 Villa Elisa, Argentina. E-mail: romero@venus.fisica.unlp.edu. ar 
Facultad de Ciencias Astronomicas y Geofisicas, UNLP, Paseo del Bosque, 1900 La Plata, Argentina. 
2 Lawrence Livermore National Laboratory, 7000 East Ave., L-413, Livermore, CA 94550, USA. 



o 
o 

Oh- 
< 
00 

> 

m 
m 
m 
m 
o 
^- 
o 



9 L,: 
6 



X 



ABSTRACT 



A recent analysis of COMPTEL data has revealed 
the existence of an unidentified and variable gamma- 
ray source, designated GRO J1411-64, at MeV ener- 
gies in the galactic plane (Zhang et al. 2002). We 
discuss the possibility that such a source can be pro- 
duced in the accreting neutron star system 2S 1417- 
624 through the bombardment of the accretion disk 
by relativistic protons accelerated in the magneto- 
sphere by the Cheng-Ruderman mechanism. 
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1. INTRODUCTION 



In a recent communication, Zhang et al. (2002) have 
reported the discovery of a new source near the galac- 
tic plane in the COMPTEL data of the Compton 
Gamma- Ray Observatory (CGRO). The source is de- 
tected at 7.2cr level in a combination of 7 viewing 
periods (VP) of the satellite, during 1995. The total 
effective exposure during Phases 1 to 5 of CGRO was 
of more than 35 days. The new source is clearly vari- 
able in the 1-3 MeV band on timescales of months. 
The maximum flux in this energy range is ~ 20 10 -5 
ph cm -2 s _1 (~ 0.3 Crab). The spectrum can be rep- 
resented by a simple power law with a photon index 
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The clear variability of the source indicates that it 
should be a compact object. Zhang et al. (2002) 
have searched for possible counterparts within the 
error location contours of the COMPTEL detection. 
There are several interesting objects in this region of 
the sky, including two high- velocity clouds (HVC), 
two radio pulsars, two Wolf-Rayet (WR) stars, one 
transient black hole binary, and a Be/X-ray transient 
source. The variability and the soft MeV spectrum 
of the source rules out the pulsars, which are typ- 
ically non- variable and hard sources, and strongly 
constrain the presumed emission region in the HVC. 



The WR stars are potential MeV emitters if they are 
colliding wind binaries. Benaglia & Romero (2003) 
discuss in detail how these kind of stars can generate 
MeV sources when they are in binary systems. In 
the present case, however, there is no evidence for 
binarity. Variable MeV emission from microquasars 
has also been recently discussed by Romero et al. 
(2002) in connection to Cygnus X-l, and remains 
as an open possibility for GRO J1411-64 (see also 
the microquasar models discussed by Romero et al. 
in these proceedings). Here, we will focus on the 
possible gamma-ray emission from accreting neutron 
star systems like 2S 1417-624, which is the Be/X-ray 
transient source within the COMPTEL error loca- 
tion contours, in an attempt to establish whether 
this is a likely counterpart. 



GAMMA-RAY PRODUCTION IN 
ACCRETING NEUTRON STARS 



Be/X-ray transients are binary systems formed by a 
massive Be star and a magnetized neutron star in an 
eccentric orbit. Be stars are rapidly rotating objects 
with strong equatorial winds. The neutron star can 
accrete directly from the circumstcllar Be wind to 
form a transient accretion disk during a periastron 
passage. Such disks have been found in systems like 
A0535+26 through the detection of quasi-periodic 
oscillations in the power spectrum of the X-ray flux 
(e.g. Finger et al. 1996a). The parameters that 
characterize the disk will change as the neutron star 
moves along the orbit. Eventually, at some distance 
that will depend on the specific system considered, 
there might be a change in the accretion regime from 
disk accretion to Bondi-Hoyle accretion. 

Once formed, the accretion disk can penetrate the 
magnetosphere of the rotating neutron star (Gosh & 
Lamb 1979). This penetration creates a broad tran- 
sition region between the unperturbed disc flow and 
the co-rotating magnetospheric flow. In this transi- 
tion zone, at the inner radius Rq the angular velocity 
is still Keplerian. Between Rq and the co-rotation ra- 
dius R co there is a thin boundary layer where the an- 
gular velocity departs from the Keplerian value. At 
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Figure 1. Sketch of the magneto sphere of an accret- 
ing magnetized neutron star. An electrostatic gap is 
produced when i7* < f2disk- Protons can be acceler- 
ated in this gap along the magnetic field lines up to 
multi-TeV energies before impacting onto the disk, 
where they produce high-energy gamma-rays through 
7T° -decays. These gamma-rays, in turn, can initiate 
electromagnetic cascades in the disk. The emerging 
gamma-rays, on the other side of the disk, can be also 
degraded in the X-ray photosphere to MeV energies 
before reaching the observer. 



R co the disc is disrupted by the magnetic pressure 
and the matter is channeled by the magnetic field 
lines to the neutron star polar cap, where it impacts 
producing hard X-ray emission. 

Cheng & Ruderman (1989, 1991) have studied the 
magnetosphere of accreting neutron stars when the 
disk rotates faster than the star, showing that three 
definite regions are formed: 1. A region coupled to 
the star by the magnetic field lines that do not pen- 
etrate the disk, 2. An equatorial region linked to the 
disk by the field attached to it, and 3. A gap en- 
tirely empty of plasma separating both regions. In 
this gap E'B / and a strong potential drop AV 
is established (see Fig. \l\ior a sketch of the situa- 
tion). For typical parameters in Be/X-ray binaries, 
AV ~ 10 13 ~ 14 V (Romero et al. 2001). Protons 
entering into the gap from the stellar co-rotating re- 
gion will be accelerated to multi-TeV energies and di- 
rected to the accretion disk by the field lines, where 
they will impact producing and ir° in pp inter- 
actions with the disk material. These pions will de- 
cay inside the disk producing electron-positron pairs, 
high-energy (TeV) photons, and neutrinos. The neu- 
trinos will escape from the system (see Anchordo- 
qui et al. 2003 for estimates of the neutrino flux), 
but the pairs and the gamma-rays initiate electro- 
magnetic cascades that will degrade the energy of 
the gamma-ray photons to form a standard cascade 
spectrum with a cut-off at MeV energies. The lep- 
tons in the disk will lose energy by inverse Compton 
(IC) interactions with the thermal X-ray photons, 



by synchrotron radiation in the field that penetrates 
the disk, and by relativistic Bremsstrahlung in the 
ions of the accreting gas. The details of the cascade 
development will be determined by the relative ra- 
tios of the cooling times among all these processes 
and the opacities for the produced photons in the 
specific cases. 



3. THE SYSTEM 2S 1417-64 

The X-ray source 2S 1417-64 was discovered by Ap- 
parao et al. (1980) using the SAS-3 satellite. Fourier 
analysis of these data revealed the presence of pul- 
sations with a period of ~ 17.63 sec. Observations 
with the Einstein satellite and optical telescopes led 
to the detection of a Be star companion (Grindlay 
et al. 1984). The distance to the system is not well 
determined due to the uncertainties in the spectral 
type; it might be in the range 1.4-11.1 kpc. The 
source was detected by BATSE (Finger et al. 1996b) 
and by RXTE (Inam et al. 2004). Large outbursts 
followed by a series of mini-bursts were reported, a 
behavior similar to what is observed in A0535+26. 
The orbital period and the eccentricity of the sys- 
tem are P = 42.19 days and e = 0.446, respectively 
(Inam et al. 2004). The X-ray data are consistent 
with the formation of a transient accretion disk dur- 
ing the outbursts. The larger outbursts last about 
two orbital periods, starting immediately after the 
periastron passage. During these outbursts the cir- 
cumstellar disk of the Be star should expand in such 
a way that it can provide matter to the neutron star 
to form a transient accretion disk when the perias- 
tron interaction occurs. The density of the circum- 
stellar disk diminishes during the mini-outbursts, in 
such a way that the accretion rate is smaller and 
the accreting material is replenished with each pas- 
sage until the disk retracts completely. The soft, 
unpulsed X-ray emission is originated in the thermal 
disk, whereas the pulsed hard X-rays are from matter 
impacting onto the neutron star polar cap. 



4. SPECIFIC MODELS 

In order to study the gamma-ray production in 2S 
1417-62 during the accretion disk phase, we adopt 
the following values for the basic parameters of the 
system: M NS = 1.4 M , R NS = 10 6 cm, B = 10 12 
G, and rj = 0.2, which correspond to the neutron 
star mass, radius, surface dipole magnetic field, and 
screening factor of the field in the disk (77 = corre- 
sponds to a perfectly diamagnetic disk), respectively. 
With these parameters and the X-ray history of the 
major outburst reported by Finger et al. (1996b) we 
calculate the structure of the accretion disk (Shakura 
& Sunyaev 1973, Gosh & Lamb 1979) and the mag- 
netosphere (Cheng & Ruderman 1991) as a function 
of time. Since the distance to the source is unknown, 
we calculated three different models corresponding 
to peak luminosities of the major outbursts of 10 , 
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Models with screening factor rj = 2 



L„„K = 1 — L p ,.k =0-5 L p „ k = 0.1 
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Figure 2. Evolution of disk parameters for three dif- 
ferent models characterized by their peak luminosi- 
ties. From left to right and from top to bottom: 1. 
Disk X-ray luminosity as a function of time (t = 
corresponds to the periastron passage), 2. Inner ra- 
dius (Ro) of the accretion disk, 3. Half-height of the 
accretion disk at Ro, 4- Particle density of the disc 
at Rq, 5. Magnetic field inside the disk at Ro, and 
6. Average energy of the disk photons at Rq. 

5 10 36 and 10 36 erg s -1 . For each model we then 
estimated the proton current impacting onto the ac- 
cretion disk surface, taking into account photo-pion 
and curvature radiation losses during the accelera- 
tion process (the proton losses resulted negligible in 
all models). Then we calculated, for different times 
of the outburst, the averaged pair cascades initiated 
in the disk and the opacities for the propagation of 
the emerging gamma-ray photons in the ambient X- 
ray fields once they escape from the disk. These re- 
sults are briefly discussed in the next section. 



5. RESULTS 

In Fig. [3 we show the main parameters that char- 
acterize the accretion disk along the evolution of a 
major outburst. The different curves in each plot 
correspond to models with different peak luminosi- 
ties, namely L pca k = 1, 0.5, and 0.1, in units of 10 37 
erg s^ 1 . In Fig. [3| we show the evolution of the 
photon density inside the disk and the energy of the 
gap-accelerated protons that impact onto the disk 
surface. These protons have energies in the ranges 
110-350 TeV and 10-70 TeV for the most and least 
luminous models, respectively. Another panel of the 
figure shows the optical depth of the disk to pro- 
ton propagation. It is clear that the disk is always 
opaque to the relativistic protons. The gamma-ray 
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Figure 3. From left to right and from top to bottom: 
1. Photon density in the disk at Rq as a function 
of time, 2. Energy of the protons impacting onto 
the disk, 3. Opacity of the disk to the protons, 4- 
Gamma-ray luminosity produced by ir° decays inside 
the disk, 5. Opacity of the disk to the ir° gamma-rays, 
and 6. The emerging ir° gamma-ray luminosity from 
the opposite side of the disk. 



luminosity produced by 7r°-decays can be as high as 
~ 5 10 35 erg s _1 for the model with i pea k = 1- 
The energy of these gamma-rays will be in the range 
20-60 TeV. The disk is mostly opaque to the prop- 
agation of these photons. The models with disks 
of smaller column densities are partially transpar- 
ent. The emerging very high energy gamma-ray flux 
(see last panel of the figure) can be detected in all 
cases by the new generation of imaging Chercnkov 
telescopes, for reasonable distances. The source, of 
course, should be observed during a major outburst. 

Those gamma-rays that are absorbed in the disk by 
photon-photon pair creation initiate electromagnetic 
cascades. Cascades are also initiated by the e^- 
pairs injected through the decay of ir^. The first 
generation of pairs mainly cools by synchrotron and 
Bremsstrahlung radiation. In Fig. 0]we show the ra- 
tio of synchrotron to Bremsstrahlung cooling times 
for these particles. Since the medium is still opaque 
to the second generation gamma-rays, there will be a 
second, and then a third generation of e ± -pairs, that 
will cool mainly by relativistic Bremsstrahlung. The 
gamma-rays that escape from the opposite side of the 
accretion disk will be injected into the X-ray photon 
field of the disk where they will be absorbed initiat- 
ing an IC cascade that will degrade their energy to 
the soft MeV band. In Fig. [5] we show the opacity 
of the X-ray field as a function of energy for t = 0, 
in the case of -L pea k = 1. This field is the addition 
of the soft blackbody disk component plus a hard X- 
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Figure 4- Ratio of synchrotron to Bremsstrahlung 
cooling times for the first generation of e^ -pairs in- 
jected into the accretion disk. 
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Figure 5. Opacity of the X-ray photosphere to 
gamma-ray propagation at t = for L poa k = 1 and 
0.1, as a function of the gamma-ray energy. 

ray component from the pulsar's polar cap that we 
have modeled as a power-law with index T = 1 and 
an exponential cut-off at kT = 200 keV. It can be 
seen that the GeV emission is completely suppressed. 
This leads to the absence of an EGRET source in this 
case. Only in the model with I/ pe ak = 0.1 the GeV 
gamma-rays can escape, and this only when the disk 
is close to its minimum luminosity. 

The MeV photons can always escape from the sys- 
tem. This emission will be variable because of the 
evolution of all relevant parameters along the or- 
bit. The final cascade spectrum will be soft, with 
an index r ~ 2.5 and a high-energy cut-off around 
~ 2m^c 4 / (£ph), where (-Eph) is the mean photon en- 
ergy of the ambient fields (see Aharonian et al. 1985 
for further details on cascades in accretion disks). 

We see, then, that several characteristics of the MeV 
gamma-ray source can be explained under the as- 
sumption that 2S 1417-64 is the counterpart: soft 
MeV spectrum, variability, and absence of EGRET 
detection. However, there are other features of the 
model that are difficult to reconcile with the ob- 
servational data: correlation with the X-ray vari- 
ability (which is expected at MeV energies, con- 
trary to what happens at GeV energies) and a ratio 
/ = Lx/L 1 > 1. In fact, for all models we expect 



/ ~ 10, instead of the observed / ~ 1. The value 
of / might be modified by a different assumption of 
the magnetic field in the neutron star, but the de- 
pendence is not strong. 



6. CONCLUSIONS 

Integral observations are crucial to establish the na- 
ture of the enigmatic source GRO J1411-64. Our 
models for the accreting neutron star 2S 1417-62 in- 
dicate that it might be a variable MeV source with a 
total gamma-ray luminosity of ~ 10 36 erg s _1 when 
the accretion disk is fully formed in models with a 
disk luminosity of 10 37 erg s -1 (a weaker source is 
predicted by the other cases studied here). Whether 
the MeV flux can be detected depends on the dis- 
tance to the source, which is unknown. If the sug- 
gested existence of anti-correlation between the MeV 
and the X-ray radiation is confirmed (Zhang et al. 
2002), 2S 1417-62 might be ruled out as counterpart 
since this anti-correlation should exist only with the 
emission above 100 MeV, due to opacity effects in the 
photosphere. In all cases a significant TeV source is 
expected. Depending on the disk density this source 
can be highly variable or roughly steady during a 
major outbursts. In some models the TeV emission 
is anti-correlated with the X-ray luminosity, as it is 
shown in the last panel of Fig. [3J Sources of this 
kind might be detectable with modern Chcrcnkov 
telescopes (Cheng et al. 1991). 
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